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Abstract 

The original proposal of Dijkgraaf, Verlinde and Verlinde for the quarter BPS dyon partition 
function in heterotic string theory on T 6 is known to correctly produce the degeneracy of dyons 
of torsion 1, i.e. dyons for which gcd(Q A P)=l. We propose a generalization of this formula 
for dyons of arbitrary torsion. Our proposal satisfies the constraints coming from S-duality 
invariance, wall crossing formula, black hole entropy and the gauge theory limit. Furthermore 
using our proposal we derive a general wall crossing formula that is valid even when both the 
decay products are non-primitive half-BPS dyons. 
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Introduction: Since the original proposal of Dijkgraaf, Verlinde and Verlinde [T] for quarter 
BPS dyon spectrum in heterotic string theory compactified on T 6 , there has been extensive 
study of dyon spectrum in a variety of Af = 4 supersymmetric string theories [21 [3j HI [HIE], 
ClElEllIDliniliailSllIlllISllISlinillB] and also in M = 8 and M = 2 supersymmetric string 
theories [T91I2T)] . However it has been realised for some time that even in heterotic string 
theory on T 6 the proposal of pQ gives the correct dyon spectrum only for a subset of dyons, - 
those with unit torsion, i.e. for which the electric and magnetic charge vectors Q and P satisfy 
gcd(Q A -P)=l [T4"ll2~T| l2"2"] . In a previous paper we proposed a general set of constraints which 
must be satisfied by the partition function of quarter BPS dyons in any N = 4 supersymmetric 
string theory and used these constraints to propose a candidate for the dyon partition function 
for torsion two dyons in heterotic string theory on T 6 [23] . In this paper we extend our analysis 
to dyons of arbitrary torsion and propose a form of the partition function of such dyons. 
Proposal for the partition function: We consider the set B of all dyons of charge vectors 
(Q, P) in heterotic string theory on T 6 , with Q being r times a primitive vector, P a primitive 
vector and Q/r and P admitting a primitive embedding in the Narain lattice [22125], he - an 
lattice vectors lying in the plane of Q and P can be expressed as integer linear combinations 
of Q/r and P. These dyons have torsion r, i.e. gcd(Q A P) = r. It was shown in [2T1I22] that 
given any pair (Q, P) of this type with the same values of Q 2 , P 2 and Q ■ P, they are related by 
T-duality transformation. We denote by d(Q, P) the index measuring the number of bosonic 
supermultiplets minus the number of fermionic supermultiplets of quarter BPS dyons carrying 
charges (Q,P) - up to a normalization this can be identified with the helicity supertrace B§ 
introduced in [26]. T-duality invariance of the theory tells us that d(Q, P) must be a function 
of the T-duality invariants, and hence has the form f(Q 2 , P 2 , Q ■ P). Then the dyon partition 
function l/<t(p, a,v) is defined as 

The sum in (TjQ) runs over all possible values of Q 2 , P 2 and Q ■ P in the set B. This in particular 
requires 

Q 2 /2 6 r 2 1, P 2 /2 e 1, Q-Per7L. (2) 

The imaginary parts of (p, a, v) in ([1]) need to be adjusted to lie in a region where the sum is 
convergent. Although the index / and hence the partition function 1/$ so defined could depend 
on the domain of the asymptotic moduli space of the theory in which we are computing the 
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partition function [13|[T8|, l22ir23j , in all known examples the dependence of $ on the domain is 
found to come through the region of the complex (p, a, v) plane in which the sum is convergent. 
Thus ([I]) computed in different domains in the asymptotic moduli space of the theory describes 
the same analytic function $ in different domains in the complex (p, a, v) plane. We shall 
assume that the same feature holds for the partition function under consideration. 

Since the quantization laws of Q 2 , P 2 and Q ■ P imply that $(p, &,v) is periodic under 
independent shifts of p, a and v by 1, 1/r 2 and 1/r respectively, eq.flT]) can be inverted as 

riMi+l/2 riM 2 +l/(2r 2 ) fiM 3 +l/(2r) 

d(Q,P) = (-l) Q - p+1 r 3 / dp da dv 



JiAh-l/2 JiM 2 -l/(2r 2 ) JiM 3 -l/{2r) 

-in(crQ 2 +pP 2 +2vQ-P) * fo\ 

provided the imaginary parts Mi, M 2 and M 3 of p, a and v are fixed in a region where the 
original sum ([T]) is convergent. 
Our proposal for a,v) is 

^(p,^)- 1 = Yl 9(s) $10 Us 2 a + ^,sv + -\ , (4) 

- -i- S k=0 1=0 ^ S S ' 



s G 2Z, s\r 



r/s 



where 



g(s) = s, (5) 

and $io(p, cr,v) is the weight 10 Igusa cusp form of Sp(2, TO). The sum over k and I in (0J 
makes $ periodic under a — > a + (1/r 2 ) and v — > v + (1/r) as required. Even though the 
function g(s) has a simple form given in ([5]), we shall carry out our analysis keeping g(s) 
arbitrary so that we can illustrate at the end how we fix the form of g(s) from the known wall 
crossing formula for decay into a pair of primitive half-BPS dyons. In particular we shall show 
that across a wall of marginal stability associated with the decay of the original quarter BPS 
dyon into a pair of half BPS states carrying charges (Qi, Pi) and (Q 2 , P 2 ) with (Qi, P x ) being 
Ni times a primitive lattice vector and (Q21P2) being N 2 times a primitive lattice vector, the 
index jumps by an amount 

Ad(Q,P)) = (-l)^*^ 

^ii|JVi J \^L 2 \N2 J 

(6) 
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for the choice g(s) = s in (jl]). Here dh(q,p) denotes the index of half BPS states carrying 
charges (q,p). When Ni = N 2 = 1 both the decay products are primitive and (EJ) reduces to 
the standard wall crossing formula [ZOEEllSllira 

Substituting (jl]) into ([3]), extending the ranges of a and v integral with the help of the sums 
over k and I, and using the periodicity of $ 10 under integer shifts of its arguments we can get 
a simpler expression for the index: 

/■iMi+l/2 piM 2 +l/(2s 2 ) pM 3 +l/(2s) 

d{Q,P) = (-l) QP+1 V(?(s)s 3 / dp da dv 

e -i« { aQZ + pP*+2vQ.P) fa s 2 & ^ ^ -1 (?) 

The set of dyons considered above contains only a subset of dyons of torsion r. This subset 
is known to be invariant under a T°(r) subgroup of the S-duality group [22J. This requires $ 
to be invariant under the transformation [23J 



a', v') = $(p, a } v) for 



p' v'\ _ (d b\ fp v\ f d c 
v 1 a 1 )~\c a)\v a) \b a ' 

a, c, d G 2Z, 6 e r ad — 6c = 1 . (8) 

On the other hand a general S-duality transformation matrix ^ a outside r°(r) will 

take us to dyons of torsion r outside the set B [22]. Thus with the help of these S-duality 
transformations on $ we can determine the partition function for other torsion r dyons lying 
outside the set E considered above. In particular if we consider the set of dyons carrying 
charges (Q', P') related to (Q, P) via an S-duality transformation 



(Q,P) = (aQ' + bP',cQ' + dP'), ( a c eSL{2, TL) 



(9) 



and denote by 1/$' the partition function of these dyons, then $' is related to $ via the relation 

p' v'\ f d b\ ( p v\ i d c 



*V,*',*') = *(A*,*) for [} r) = {~ a) a)' W 

This allows us to determine the partition function of all other sets of torsion r dyons from 
the partition function given in In particular for r = 2, choosing ^ = ^q l^) we 
recover the dyon partition function proposed in 
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We shall now show that the proposed partition function (TJJ satisfies various consistency 
tests described in [23] . 

S-duality invariance: We shall first verify the required S-duality invariance of the partition 
function described in eq.flU). Using Sp(2, 7L) invariance of $10 (x, y, z), and that b is a multiple 

of r for f a , 1 G T°(r) one can show that 



c d 



* UI ( p', sV + 4, sv' + - ) = $10 ( P, s 2 cr + ^, st; + - ) . ( I I ; 



where 

k' = kd 2 - 2cdlr G "Z, V — (ad + bc)l — bdk/r G 2£ . (12) 

Thus 

E E sV + 4. + i) = E E ^ U s2& + S> s?} + > ( 13 ) 

fc=0 i=0 ^ ' fe'=0 Z'=0 ^ ' 

and we have the required relation (jSJ). 

Wall crossing formula: We shall now verify that (j7|) is consistent with the wall crossing 
formula. As in [23] we shall only consider the decay into a pair of half-BPS dyons [T3"ll3"5'| l3"(3 | l3"T] 

(Q,P)^(Q 1 ,P 1 ) + (Q 2 ,P 2 ), (14) 

(Q x , Px) = (aQ + /3P, 7 Q + 5P), (Q 2 , P 2 ) = (5Q - (3P, -jQ + aP) , (15) 

aS = (3j, a + 5 = l. (16) 

Since any lattice vector lying in the plane of Q and P can be expressed as a linear combination 
of Q/r and P with integer coefficients, we must have /3, 5, a G 2, 7 G 22/r. Thus we can 
write 7 = j'/K, where K £ IK, K\r and gcd(7', iC)=l. The condition a;5 = /? 7 together with 
the integrality of at, (3,5 now tells us that (3 must be of the form K(3' with (3' G Z. Thus we 
have 

= 7 = ^, K,a,8,p,ie 2Z, tf|r, gcd( 7 ',^) = l. (17) 

Using eqs. flIB"]) . (TTTT) we have 

a + 5 = 1, a6 = Pi, a : (3',^,5e TL. (18) 

The analysis of [13J now shows that we can find a', b', c', d' such that 

a = a'd', (3' = —a'b', 7' = c'd', 5 = —b'c', a' , b', c', d' G ad' — b'c' = 1 . (19) 
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As a consequence of (fT9l) and the relation gcd(7', K) = 1 we have 

gcd(a', b') = gcd(a', c') = gcd(c, d!) = gcd(6', d!) = 1, gcd(c', K) = gcd(d' } K) = 1 . (20) 
Using eqs. (jIT|) - (TT9~]) we can now express (jT5l) as 

(Qi.Pi) = (a'K,c')(d'KQ/r-b'P), (Q 2 ,P 2 ) = (b'K,d')(-c'KQ/r + a'P), K = r/K . 

(21) 

Since according to (TSUI) . gcd(a'i^, c')=l, gcd(6'if, c?')=l, and since any lattice vector lying in 
the Q-P plane can be expressed as integer linear combinations of Q/r and P, it follows from 
(I2TI) that (Qi, Pi) can be regarded as Ni times a primitive vector and (Q 2 , P2) can be regarded 
as N 2 times a primitive vector where 

N\ = gcd(d'K, b') = gcd(K, b'), N 2 = gcd(c'K, a) = gcd(K, a) . (22) 

In the last steps we have again made use of ( 1201) . It follows from ( 1201) and ( 1221) that 

gcd(N 1 ,N 2 ) = l, N X N 2 \K. (23) 

We shall now use the formula (J7J) for the index in different regions of the moduli space and 
calculate the change in the index as we cross the wall of marginal stability associated with the 
decay (TH1) . For this we need to know how to choose the imaginary parts of p, a and v along 
the integration contour for the two domains lying on the two sides of this wall of marginal 
stability. A prescription for choosing this contour was postulated in [23] according to which as 
we cross the wall of marginal stability associated with the decay (T141) . the integration contour 
crosses a pole of the partition function at 

p~f-a(3 + v(a-5) =0. (24) 

Thus the change in the index can be calculated by evaluating the residue of the partition 
function at this pole. We shall now examine for which values of s the $io(p, s 2 ®-, sv)~ l term in 
the expansion ([7]) has a pole at (l24l) . The poles of $io(p, s 2 a, sv) -1 are known to be at 

n 2 s 2 (pa — v 2 ) + riis 2 a — mip + m 2 + jsv = 

j 2 1 

m 1 ,ni,m 2 ,n 2 e TL, j e2 2K + 1, m\n\ + m 2 n 2 + — = - . (25) 
Comparing (1251) and (1241) we see that we must have 

m 2 = n 2 = 0, j = -(a-5), n 1 = -^-P, mi = -\j , (26) 

s s z 
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for some A. The last condition in (J25|) . together with eqs. tjlb]) now gives 

A = a . (27) 

Thus we have from ffTTj) . (|26|) 

j = a - S, m l = -s 7 = -is/K, n x = -j3/s = -Kf3'/s . (28) 
Since gcd(7', K) = 1, the second equation in fl28|) shows that s must be a multiple of i^: 

s = K§, sell. (29) 

Substituting this into the last equation in (1281) and using (|T9|) we see that 

a'V a'b' m , . 

m = G TL. (30) 

s s 

Since gcd(a', b')=l, we must have a unique decomposition 

a = LiL 2 , Li\b', L 2 \a' . (31) 
On the other hand since s divides r, it follows from (!29l) that s must divide r/K — K. Thus 

L^K, L 2 \K. (32) 



It now follows from (1221) that 

Lx\N x , L 2 \N 2 . (33) 

Conversely, given any pair (L 1; L 2 ) satisfying (|33l) . it follows from (|22|) that L x , L 2 will satisfy 
(EH), ( |32l) . This allows us to find integers mi, ni, j satisfying ( 1261) via eqs.(|271 ) -( l3~ll) . 

This shows that the poles of $(p, a, w) _1 at ff24|) can come from the s = KL\L 2 terms in 
(jlj) for Li|iVi and L 2 \N 2 . Our next task is to find the residues at these poles to compute the 
change in the index as we cross this wall. For this we define 

a = a'/L 2 , d = d'L 2 , b = b'/L u c = c'L u s = KL X L 2 . (34) 

It follows from ( 13~TT) that a, b, c, d are all integers. In terms of these variables the location of 
the pole given in ( 1241 can be expressed as 

Sq 1 cdp + dbs a + {ad + bc)v = . (35) 
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We now define 



p = d 2 p + b 2 sl& + 2bds v, a' = c 2 p + a 2 sl& + 2acs v, v' = cdp + abs^a + (ad + bc)s v . (36) 

The change of variables from (p, Sq<7, sqv) to (p',a',v') is an Sp(2, TL) transformation. Thus 
we have 

$io(p, s 2 a, s v) = $i (p, a', v') . (37) 
In the primed variables the desired pole at ( 1351) is at v' = 0. We also have 

pP 2 + dQ 2 + 2vQ-P = p'P' 2 + a'Q' 2 + 2vQ' ■ P' , (38) 

where 

Q' = d Q/s - bP, P' = -cQ/sq + dP . (39) 

Finally we have 

dp'da'dv' = sldpdvda . (40) 

This is consistent with the fact that <3>io(p', <r', v') is invariant under integer shifts in p', a' 
and v' so that in the primed variables the volume of the unit cell is 1, while in the unprimed 
variables the volume of the unit cell is 1/sq. We can now express the change in the index from 
the s = So term in (J7J) as 

(Ad(Q,P)) ao = (-lf- p+1 g(s ) / dp' ' da'idv' 



JiM[-l/2 JiM^-1/2 J 

e -™(*'Q' 2 +p'p' 2 +wQ'-n $ 10 & ' : v'y 1 , (41) 

where the v' contour is around the origin, - as in [23] we shall use the convention that the 
contour is in the clockwise direction. Using the fact that 

d> 10 (p\ a', v') = -47T 2 v' 2 V (p') 2A V(*T + O (v' A ) , (42) 

near v' = 0, we get 

«iM[+l/2 



(Ad(Q,P)) S0 = (-l) Q - p+1 g(s )Q'-P' dp'e-^' p ' 2 V (pT 24 

JiM[~l/2 

riM^+l/2 

/ da'e-^' QI \{a')- U 

J iM' -1/2 



iM' 2 -l/2 

(-1) Q P+1 g(s ) Q' ■ P' d h (Q', 0) d h (P', 0) , (43) 
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where dh(q, 0) denotes the index measuring the number of bosonic half BPS supermultiplets 
minus the number of fermionic half BPS supermultiplets carrying charge (q, 0). 

We shall now express the right hand side of ( 1431) in terms of the charges (Q\,P\) and 
{Q2, P2) of the decay products. First of all it is easy to see that 

- (_i)Qi- p 2-Q 2 --Pi _ 

Furthermore it follows from (12TT) . (|34|) and (1391) that 

(Qi, Pi) = L x (a'KQ', c'Q') , (Q2, P2) = L 2 (b'KP', d'P') , (45) 

Qi-P2-Q2-Pi = s Q'-P'. (46) 

Now according to (1201) the pair of integers (a'fT, c') are relatively prime, and the pair of integers 
(b'K, d') are also relatively prime. Thus using S-duality invariance we can write 

4 (f?Zr) = 4(aW,cV)=4(<2',0), 

4 (§'£) = d hiP'KP\d'P') = d h {P'^). (47) 
Using dHJ, (Jl6l) and (j4"T|) we can now express (j4"3j) as 



(Ad(Q,P)) S0 = (-l)Q^-Q^g( So )-(Q 1 -P 2 -Q 2 -P 1 )d h (%^) d h 

so V-^i L iJ 



Q2 P2 

L 2 L 2 



(4£ 



For a given decay K is fixed but Li and L 2 can vary over all the factors of N\ and N 2 . Thus 
the total change in the index is obtained by summing over all possible values of so of the form 
KL\L 2 . Thus gives 

Ad(Q,P)) = (-i)^ p ^ p ^(Q 1 .p 2 -Q 2 .p l ) 

EE^)^4(||)4(||). (49, 

We can now fix the form of the function g(s) by considering a decay where the decay products 
are primitive, i.e. Ni — N 2 — 1. In this case we have L\ — L 2 — 1, and ()49l) takes the foruQ 

Ad{Q, P)) = (Qi . p 2 _ g 2 . Pl) ^ _L 4 ( Ql) Pl) 4 (g 2; p 2 ) . (50 ) 



1 In this argument we have implicitly assumed that for a given K, it is possible to find integers a' , b' , c' , 
d' satisfying a'd' - b'c' = 1, gcd(c',if) = gcd(d',K) = gcd(a',K) = gcd{b',K) = 1, so that ^ holds and 
we have Ni = N2 — 1 according to (f22|) . If cither K or K is odd then this assumption holds with the choice 
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In order that this agrees with the standard wall crossing formula for primitive decay 
ESEQiEUdniESIGS] we must have g(K) = K. Since this result should hold for all K\r we see 
that we must set g(s) = s as given in (jBJ). Using this we can simplify the wall crossing formula 
fj49l) for generic non-primitive decays to the form given in ([6]). 

Black hole entropy: In order to reproduce the leading contribution to the black hole entropy 
in the limit of large charges, the partition function must have a pole at [T] 

pa - v 2 + v = 0. (51) 

Furthermore, in order to reproduce the black hole entropy to first non-leading order, the inverse 
of the partition function near this pole must behave as 



$(p, cr, v) oc (2v - p - a) 10 {v 2 V (p) 24 r/(a) 24 + 0(v 4 )} , (52) 

where 

pa - v 2 pa-(v-l) 2 pa-v 2 + v . . 

p= : , o = : , v = : . (53) 

a a a 

We shall now examine the poles of $ given in (j3J) to see if it satisfies the above relations. For 

this we recall eq. fl25|) giving the pole of $io(/5, s 2 a, si))' 1 . Comparing fl25|) with ([ST]) we see 

that in order to get a pole at flBT]) we must choose in (12"5"1) 

n 2 = \/s 2 , j = X/s, n\ = mi = m 2 = , (54) 

for some A. The requirement m\ri\ + m 2 n 2 + \j 2 = \ then gives 

A = s, 7i2 = - ■ (55) 
s 

Since n 2 must be an integer this gives s = 1. Thus only the s = 1 term in (jlj) has a pole at 
(IBTI) . It follows from the known behaviour of $iq near its zeroes that $ defined in (j3J) satisfies 
the requirement (1521 . 

Gauge theory limit: Finally we shall show that in special regions of the Narain moduli space 
where the low lying states in string theory describe a non-abelian gauge theory, the proposed 



If K and K are both even then we cannot find a', b' , c' and d! satisfying 



all the requirements since in order to satisfy a'd! — b'd = 1 at least one of them must be even. However in this 

case if we choose = \ \ \ \ ^ en we satisfy (|20| and have N\ — 1 and N2 — 2 according to (|22|) . 

We can now demand that the wall crossing formula in this case agrees with the one derived in 23J for decays 
where one of the decay products is twice a primitive vector. This gives g{K ) — K even for K, K both even. 
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dyon spectrum of string theory reproduces the known results in gauge theory. Since the T- 
duality invariant metric L in the Narain moduli space descends to the negative of the Cartan 
metric in gauge theories, and since the Cartan metric is positive definite, all the gauge theory 
dyons have the property 

Q 2 <0, P 2 <0, Q 2 P 2 > (Q ■ P) 2 . (56) 

Thus in order to identify dyons which could be interpreted as gauge theory dyons in an appro- 
priate limit we must focus on charge vectors satisfying (1561) . 

In order to identify such dyons we need to expand the partition function $ _1 in powers of 
e 2m P, e 2mcr and e 2mv and pick up the appropriate terms in the expansion. For this we need to 
identify a region in the (p, a, v) space (or equivalently in the asymptotic moduli space) where 
we carry out the expansion, since in different regions we have different expansion. We shall 
consider the region 

3(p), 3(a) » -%v) » 1 , (57) 

where 3(z) denotes the imaginary part of z. The results in other relevant regions will be 
related to the ones in this region by S-duality transformation. In the region (1571) the only term 
in $io (p, s 2 a + ^,sv+ |) 1 which has powers of e 27n/5 , e 2ntcr and e 2mv compatible with the 
requirement (1511 is 

e -2nip-2Tri(s 2 &+-$ I )-2ni(sv + t) 

This is in fact the leading term in the expansion in the limit (}5Tj) . Substituting this into (J3|) 
and performing the sum over k and I we see that only the s = r term in the sum survives and 
the result is 

r e ~2iri(p+r 2 a+rv) 

This corresponds to dyons with 

Q 2 /2 = -r 2 , P 2 /2 = -l, Q-P = -r, (60) 

with an index of (— l) r+1 r. 

We can also determine the walls of marginal stability which border the domain in which 
these dyons exist. This requires determining the region in the (S(p), S(<5"), Q^v)) space in which 
the expansion (|59|) of $io(p, r 2 a, rv)" 1 is valid since then we can determine the associated walls 
of marginal stability using (12~21) . For this we shall utilize the known results for r = 1; in this 
case the walls of marginal stability bordering the domain in which fl59l) is valid correspond to 
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the decays into (Q, 0) + (0, P), (Q - P, 0) + (P, P) and (0, P - Q) + (Q, Q) respectively [38]. 
Using ( 124"I) we now see that validity of the expansion of $io(/5, v)' 1 given by ( 1591) with r = 1 
is bounded by the following surfaces in the (3(/5), 9(v)) space: 

%(v)=0, %{v + &) = 0, S(v + p) = 0. (61) 

We can now simply scale a by r 2 and # by r to determine the region of validity of the expansion 
fl59l) for $i (p, r 2 a,rv): 

$(v) = 0, Q(v + ra) = 0, %{rv + p) = . (62) 

Comparing these with f[2"4"|) we now see that the corresponding walls of marginal stability are 
associated with the decays into 

(Q, 0) + (0, P), (Q- rP, 0) + (rP, P), (o, P - V) + K?, ^Qj . (63) 

Let us now compare these results with dyons in M = 4 supersymmetric SU (3) gauge theory. 
If we denote by «i and a pair of simple roots of SU(3) with a\ = a 2 . = —2 and ai ■ a2 = 1, 
then the analysis of [39lll0llll],ll2],|l3] shows that the gauge theory contains dyons of charge 

(Q,P) = K-oj), (64) 

with index (— l) r+1 r. These are precisely the dyons of the type given in fl60l) . Furthermore 
using string junction picture [HISS], ref. [3H] also determined the walls of marginal stability 
bordering the domain in which these dyons exist. These also coincide with fl6"31 . 

The spectrum in gauge theory contains other dyons of torsion r related to the ones described 
above by S-duality transformation. Since our construction is manifestly S-duality invariant, 
the results for these dyons can also be reproduced from the general formula given in (T41). 

Gauge theory also contains other dyons which are not related to the ones described here by 
S-duality [4X114*2] . These typically require higher gauge groups and has additional fermionic zero 
modes besides the ones required by broken supersymmetry. Quantization of these additional 
zero modes gives rise to additional bose-fermi degeneracy, and as a result the index being 
computed here vanishes for these dyons. This is also apparent from the fact that these dyons 
typically exist only in a subspace of the full moduli space; as a result when we move away from 
this subspace the various states combine and become non-BPS. Some aspects of these dyons 
have been discussed recently in [4B1I43]. 
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